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Abstract 



The properties of high-density nuclear and neutron matter are studied using 
a relativistic mean-field approximation to the nuclear matter energy func- 
tional. Based on ideas of effective field theory, nonlinear interactions between 
the fields are introduced to parametrize the density dependence of the en- 
ergy functional. Various types of nonlinearities involving scalar-isoscalar (a), 
vector-isoscalar (u), and vector-isovector (p) fields are studied. After calibrat- 
ing the model parameters at equilibrium nuclear matter density, the model 
and parameter dependence of the resulting equation of state is examined in 
the neutron-rich and high-density regime. It is possible to build different mod- 
els that reproduce the same observed properties at normal nuclear densities, 
but which yield maximum neutron star masses that differ by more than one 
solar mass. Implications for the existence of kaon condensates or quark cores 
in neutron stars are discussed. 
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I. INTRODUCTION 



Numerous calculations have established that relativistic mean-field models provide a 
realistic description of the bulk properties of finite nuclei and nuclear matter In 
addition to this successful low-energy phenomenology, these models are often extrapolated 
into regimes of high density and temperature to extract the nuclear equation of state (EOS), 
which is the basic ingredient in many astrophysical applications and in microscopic models 
of energetic nucleus-nucleus collisions. 

Based on the original version of Walecka [[| and its extensions [|]||, relativistic mean- 
field models generally involve the interaction of Dirac nucleons with neutral scalar and 
vector mesons and with isovector p mesons. One of the key observations in their success 
is that to provide sufficient flexibility, nonlinear self-interactions for the scalar meson must 
be included P,[^J^-|IT|. Since these models were proposed to be renormalizable, the scalar 
self-interactions are limited to a quartic polynomial, and scalar-vector and vector-vector 
interactions are not allowed ||12| . One of the motivations for renormalizability, as discussed 
in Walecka's seminal paper, is that once the model parameters are calibrated to observed 
nuclear properties, one can extrapolate into regimes of high density or temperature without 
the appearance of new, unknown parameters. 

An alternative approach is inspired by effective field theories, such as chiral perturba- 
tion theory Hl3| , |14|1 , which successfully describes the low-energy phenomenology of hadronic 



Goldstone bosons |I7)| , [TE|] . Although a lagrangian usually serves as the starting point, the 
meson and baryon fields are no longer considered elementary, and the constraint of renor- 
malizability is dropped. This has several important consequences. First, there is no reason 
to restrict meson self-interactions to a simple quartic polynomial in the scalar field; on the 
contrary, one should include all interaction terms that are consistent with the underlying 
symmetries of QCD. Second, since there are an infinite number of coupling constants, one 
must find suitable expansion parameters for the systems under consideration, and one must 
develop a systematic truncation scheme for the effective theory to have any predictive power. 
Third, extrapolation of calculated results into new regimes of the physical parameters be- 
comes problematic, because the truncation scheme may break down, and predictions can 
become sensitive to unknown parameters. 

Within the framework of effective field theory, mean-field models of nuclear structure 
and the EOS must be interpreted in a new context. One important observation is that near 
normal nuclear density, the mean scalar and vector fields (or nucleon self-energies) , which we 
denote as $ and W, are large on nuclear energy scales but are small compared to the nucleon 
mass M and vary slowly in finite nuclei. This implies that the ratios $/M and W/M and 
the gradients |V$|/M 2 and IVW^I/M 2 are useful expansion parameters. The assumption of 
"naturalness" in effective field theory is also important. Naturalness implies that the coef- 
ficients of the various terms in the lagrangian, when expressed in appropriate dimensionless 
form, should all be of order unity. When combined with meaningful expansion parameters, 
this means that one can anticipate the approximate magnitude of mean-field contributions 
to the energy (at least up to moderate nuclear densities) and thereby motivate a suitable 
truncation scheme; if the coefficients are natural, the omitted terms will be numerically 
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unimportant .[] Naturalness also implies that one should include all possible terms (that is, 
those allowed by the symmetries) through a given order of truncation; it is unnatural for 
some coefficients to vanish without a relevant symmetry argument. 

From this point of view, it is difficult to justify nuclear mean-field models that include 
only scalar self- interactions |4||||7],|1(J , and recently, generalizations that also include quartic 
self- interactions for the neutral vector meson have been discussed P, |T9l , |20|| . Moreover, a 



new analysis involving all meson self-interactions through fourth order in the isoscalar scalar 
and vector fields has been performed ]2l|. These extensions give rise to additional model 



parameters (coupling constants) that must be constrained by calibrating to observed nuclear 
properties. For the truncation at fourth order to be sensible, the parameters so obtained 
should exhibit naturalness. 

Although it is possible to discuss effective hadronic field theory from the point of view of 
a lagrangian, as above, the expansion in powers of the mean fields is a low-density expansion, 
and it is hard to justify the neglect of many-body corrections, which are known to be relevant 
in nuclear structure and in the EOS. Alternatively, one can consider this expansion at the 
level of an energy functional or effective action [^TJ. In such a formulation of the relativistic 
nuclear many-body problem, the central object is an energy functional of scalar and vector 
densities (or more generally, vector four-currents) |22|-21|. Extremization of the functional 



gives rise to Dirac equations for occupied orbitals with local scalar and vector potentials, not 
only in the Hartree approximation, but in the general case as well. Rather than work solely 
with the densities, one can introduce auxiliary variables corresponding to the local potentials, 
so that the functional depends also on mean meson fields. The resulting field equations have 
the same form as in a Dirac-Hartree calculation [[H]], but correlation effects can be included, 
if the proper energy functional can be found. This procedure is analogous to the well-known 
Kohn-Sham JZH] approach in density-functional theory, with the local meson fields playing 
the role of (relativistic) Kohn-Sham potentials; by introducing nonlinear couplings between 
these fields, one can implicitly include additional density dependence in the potentials. Thus 
the nonlinear meson interaction terms simulate more complicated physics, such as one- and 
two-pion exchange, or vacuum-loop corrections, which might be calculated directly in a 
more microscopic many-body approach p0| , |2^] . The fields (and their gradients) again serve 
as useful expansion parameters at moderate density, so the nonlinear interaction terms can 
be truncated, leaving a finite number of unknown couplings. 

Rather than focus on the calculation of the nonlinear couplings from an underlying 
effective lagrangian, we wish to concentrate instead on how well the energy functional can 
be calibrated by fitting the couplings to observed nuclear properties, and on the limitations 
on the extrapolation of the resulting EOS into the high-density regime. In general, even 
with a significant truncation, the number of unknown couplings exceeds the number of 
normalization conditions, which we take to be five properties of infinite nuclear matter: 
the equilibrium density and binding energy (p , — e ), the nucleon effective (or Dirac) mass 
at equilibrium (Mq), the compression modulus (K ), and the bulk symmetry energy (a 4 ). 
(Experience has shown that an accurate reproduction of these five properties leads to realistic 



x It has also been shown recently that the naturalness assumption is consistent with dimensional 
counting in chiral perturbation theory JlTjJl^] . 
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results when the calculations are extended to finite nuclei IPJ7I J10IJTTIJ2TII .) Thus families of 
models can be generated which describe exactly the same nuclear matter properties at 
equilibrium [HJ. We then investigate the differences in the high-density EOS predicted by 
models within a given family. This question is important for astrophysical applications, 
particularly in the study of neutron stars; it will be difficult to deduce the existence of 
"exotic" neutron-star structure (for example, hyperonic matter, kaon condensates, or quark 
cores) unless the EOS of the more mundane components (neutrons, protons, and electrons) 
is well constrained. 

If our truncation of the energy functional is motivated by low-density behavior, why 
should we have any confidence at all in a high-density extrapolation? This is indeed the 
crucial question, and we are not attempting here to justify such an extrapolation; we are 
merely recognizing that this procedure is often used in neutron-star calculations, even recent 
ones, without any mention of the implicit assumptions about the absence of additional con- 
tributions at high density pj| -p!^ . We therefore feel it is timely to investigate quantitatively 
the uncertainties in the extrapolated equation of state. 

We begin the theoretical analysis with a model that contains meson-meson and meson 
self-interactions described by an arbitrary finite polynomial in the fields. We find that both 
the asymptotic (high-density) limit of the EOS and the approach to this limit (i.e., the 
"stiffness") are model dependent. In particular, one can construct models with the same 
equilibrium nuclear matter properties that yield high-density equations of state ranging all 
the way from the causal limit (p = £) to one that resembles a free relativistic gas (p = £/3) 
PJl9|. (Here p is the pressure and £ is the energy density.) 

As an explicit example, we consider a model that includes self-interactions for the 
isoscalar scalar and vector mesons and for the p meson up to fourth order in the fields. 
To our knowledge, this is the first time nonlinear terms in the p meson mean field have 
been included. (Note that the p field enters here as an effective field whose purpose is to 
parametrize the isospin dependence; thus, fundamental questions about causal propagation 
|33| , |34|1 and spin mixing []35] are not relevant.) To provide a quantitative measure of the 



variations in the EOS, we compute neutron star masses, which turn out to be sensitive to 
the model types and to changes in the parametrizations, even for models that reproduce 
the same equilibrium nuclear matter properties. In some cases, variations in the calculated 
maximum mass are more than one solar mass. 

These results lead us to two basic conclusions. First, existing methods for calibration of 
the EOS at normal density are not sufficient to provide a satisfactory extrapolation into the 
density regime relevant for neutron stars. The basic problem rests with the quartic neutral 
vector meson (u) interaction, which produces major modifications in the high-density behav- 
ior. We discuss the existing situation regarding the specification of this term and prospects 
for improved calibration. Second, we find that if the quartic u term can be accurately de- 
termined, additional higher-order terms (and the quartic p term) produce relatively minor 
changes in the neutron star mass. This occurs because the W 4 term softens the high-density 
EOS so completely that additional interactions have little effect, at least in the density 
regime relevant for neutron stars. 

These conclusions have important implications, because the contributions to the EOS 
from the neutrons and protons are commonly believed to be the best understood part of the 
physical input necessary to describe a dense stellar object. Without well-constrained results 
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from this part of the EOS, it will be impossible to determine the importance of additional 
dynamics, for example, the transition from nuclear matter to quark matter in the interiors 
of neutron stars p6|j37|j27|] , or the role of strangeness in the form of hyperons p8| , |39| or a 
kaon condensate 
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As an example of how uncertainties in the basic nuclear EOS can 
influence these interesting effects, we study the transition from hadronic matter to quark 
matter using a simple model |37|,p| and show the possible variations in the results. 

The outline of this paper is as follows: In Sec. II, we present the general model, which 
involves an arbitrary number of nonlinear meson interactions, and derive the EOS. Based on 
this general model, Section III is devoted to the high-density limit of the EOS. In Sec. IV, 
we apply our model to neutron stars. For the quantitative analysis, we initially include self- 
interactions up to fourth order in the p and u fields and then investigate the consequences 
of sixth-order and eighth-order uo self-interactions. We also briefly discuss the parameter 
and model dependence of the transition to quark matter in the central region of the star. 
Section V contains a short summary and our conclusions. 



II. THE NUCLEAR EQUATION OF STATE 

We describe the nuclear equation of state using a relativistic approach involving valence 
Dirac nucleons and effective mesonic degrees of freedom, which are taken to be neutral 
scalar and vector fields, plus the isovector p meson field. Rather than focus directly on a 
lagrangian, we consider instead as a starting point an effective action 

r = r[0,^,b M ] , (i) 

which is a functional of the meson fields denoted by 0, V^, and b M for the scalar, vector- 
isoscalar, and vector-isovector field, respectively. In principle, this functional can be calcu- 
lated in a many-body approach based on a Lagrangian for the nucleon-nucleon interaction, 
or its general form might be obtained from an underlying theory. Here we will be satisfied 
to parametrize the effective action, calibrate it as accurately as we can to observed nuclear 
properties, and then examine the predicted high-density equation of state. 
The effective action is related to the thermodynamic potential Q |4]J by 



^ = r[0,^,b M ] , (2) 

where the fields are determined by the general thermodynamic principle that they should 
make Q stationary: 



A basic property of the functional is that it reflects the underlying symmetries [f£q| . Thus, if 
we assume that the system possesses two conserved charges, namely, baryon number B and 
the third component of total isospin I 3 , this gives rise to two chemical potentials /i and v: 

n = n(p,fji,u), (4) 

with 
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h = -J d ^ Ps= --, (6) 

/? the inverse temperature, and p 3 = p p — p n . Note that the fields can be held fixed in 
evaluating the partial derivatives in Eqs. ([5]) and (||) by virtue of the extremization conditions 



In the field theoretical context [43|, one can show that the effective action can be ex- 
panded as a power series in gradients of the fields. Thus, after taking the zero-temperature 
limit, one usually writes 

r[0,V„,b„] = J rf 4 x[-V cfr (0,^,b M ) + Z(<f>,V /t ,b ft ,d v (f>,d v V li ,d v b /t )] , (7) 

where the second term vanishes in a uniform system, and the dependence on the chemical 
potentials has been suppressed. For the effective potential, we make the following Ansatz: 

V eff (0, V„ b M ; W) = \ ™ s 2 2 - \ m^/" - \ mJb M -b^ 

+ itrlnZ^,(/i, v) — itrlnZ^O, 0) 
+ AV(0,V^,b,-b>,z/) , (8) 

with a nonlinear potential 

av(0, v^w, k-W; ^u) = -J2 <k jh {ji, v) <p (v^Y (K-bn k (9) 

i,j,k 

that contains at least three powers of the fields: i + 2j + 2k > 3. The fermionic contributions 
are represented by a one-body term tr In (with the appropriate zero-density subtraction 
f44| ) and by terms where the fermions have been "integrated out", which results in a (gen- 



erally nonanalytic) dependence of the mesonic coefficients on the chemical potentials /i 
and v. There will also be contributions to the that are independent of the chemical po- 
tentials; these arise from integrating out heavy degrees of freedom and vacuum loops. Thus 
V e ff contains explicit contributions only from valence nucleons and classical meson fields. 
The fermionic part tr In is obtained by evaluating the trace of the kernel 

Kiji, u) = (id" - g v V» - l - g p rW) llx + /i 7 ° + \^ — (M — g s 0) (10) 

using Dirac wave functions calculated in the presence of static background fields. The 
subtraction removes contributions from negative-energy states, which are already included 
implicitly in the nonlinear parameters a^. (See Ref. [13] f° r an analogous calculation.) The 



values of the fields are determined by extremization at the given values of [i and v. 

The potential of Eq. (|8p represents an effective field theory for the interacting nucleons. 
Although the exact form of the effective potential is not known, we introduce the meson 
mean fields as relativistic Kohn-Sham potentials and consider the valence nucleons 
moving in the resulting local fields. The nonlinear interactions of the fields generate implicit 
density dependence above and beyond that arising from the couplings in Eq. (|T0|) . Thus the 
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series in Eq. (Qj can be interpreted as a Taylor series parametrization of the unknown part 
of the effective potential, which includes the effects of nucleon exchange and correlations, as 
well as contributions from other mesons and the quantum vacuum. (See the discussion in 



Refs. 21.44 



Although the couplings generally depend on the chemical potentials, experience with 
calculations for finite nuclei and nuclear matter, together with explicit computations of ex- 
change and correlation corrections |]45|| , implies that mean fields and constant couplings 
provide an adequate (albeit approximate) parametrization of these many-body effects. Thus 
we consider the as constants in the sequel and leave the study of their dependence on fi 
and v as a topic for future investigation. Moreover, at low densities and temperatures, the 
mean meson fields are small compared to the nucleon mass, and so provide useful expansion 



parameters |21|] . Thus, in practice, the series in Eq. ([)]) can be truncated at some reasonable 
order, and the relevant question in this paper is how far one can extrapolate the truncated 
potential into the high-density regime. 

In principle, the unknown coefficients (coupling constants) can be constrained by im- 
posing chiral symmetry and other symmetries of the underlying QCD, such as broken scale 
invariance. (Lorentz covariance and isospin symmetry are already incorporated explicitly.^) 
As has been discussed recently, however |TT|f2T|], if one assumes a nonlinear realization of the 



chiral symmetry for the pions and nucleons > the meson interaction terms are essentially 
unconstrained.^ We therefore take the couplings as free model parameters in our approach. 
Obviously, an infinite number of normalization conditions is generally needed to fix their 
values. Since this is not feasible in practice, we terminate the summation at the finite values 

^maxj jmaxj and k max . 

According to Eq. ([3D, the thermodynamic potential Q must be stationary with respect 
to changes in the fields for fixed values of the proton and neutron chemical potentials 

^ = A* + \ = (4 P + M* 2 ) 1/2 + W + \R , (11) 
»n = » - \ = (4„ + M* 2 ) 1/2 + W--R. (12) 



Here, following Bodmer fL9|| , we define the scaled meson fields $ = g s (p, W = g v Vo, and 
R = g p bo, with bo the timelike, neutral part of the p meson field; the effective nucleon mass 
is M* = M — $. (We work in the rest frame of the infinite matter, where the spatial parts 
of the vector fields vanish.) The Fermi momenta for protons (fcpp) and neutrons (k^ a ) are 
related to the conserved baryon density 

1 

37 2 



P=^( fc Fp + fcFn) (13) 



2 Note that since the energy functional is an effective functional, we presently know of no reason to 
exclude terms that explicitly contain the medium four-velocity u^, such as u^V^VVv This issue 
will be considered in a later publication. 

3 Broken scale invariance leads to restrictions on the purely scalar interactions, as shown in 
Refs. [47 1 and p3[], but we will not consider these limitations here. As discussed in the next 



section, the details of the scalar dynamics do not play a major role in our analysis. 



7 



and isovector density 



P*= 3^{ k h- k Fn)- ( 14 ) 



Using Eqs. (Q), (0), and (^), together with relations (|TTJ) and ( |12"D for the chemical 
potentials, it is straightforward to eliminate the chemical potentials in favor of the densities 
and to compute the pressure p and the energy density £ = —p + \ip + \up- A : 

1 /- fc F P k 4 1 f**n „ 



3vr 2 7o (Jfe 2 + M* 2 ) 1/2 3tt 2 Jo (P + M* 2 ) 1 ^ 
+ -^^ + -Li? 2 --L$ 2 + E^^^^, (15) 

ZC v zc p ZC S jj,fe 

£ = i_ / fcFp dfc fc 2 (fc 2 + M* 2 ) 1 ' 2 + -i / feFn dA; A; 2 (A; 2 + M* 2 ) 1 ' 2 

7T 2 JO 7T 2 JO 

+ Wp+ I^ 3 _ 1 1 1 $ 2_^^.. fc$ i K ,2 i ^ (16) 

v p s 

Here the ratios c 2 = g 2 /m 2 and = a^k/ {g\g 2 ^ g 2k ) have been introduced for convenience. 

As noted earlier, the pressure in Eq. flI5]) and the energy density in Eq. ( |i6|) also contain 
vacuum contributions arising from the partition function of the nucleons. However, at least 
at the one-baryon-loop level, these vacuum terms can be absorbed in the definition of the 



nonlinear couplings in Eq. (gj) f44jl ) and thus we include explicitly only the contributions 
from valence nucleons. 

At zero temperature, the stationarity conditions of Eq. (|3]) with fixed chemical potentials 
are equivalent to an extremization of the energy at fixed baryon and isovector density. This 
leads to the self-consistency equations 

\ $ - £ z a ijk W 2 > R 2k = p s , (17) 

° s i,j,k 

W + ]T 2j a ijk & W 2 '- 1 R 2k = p , (18) 
I + £ 2A; a ijk & W 2j R 2k ~ l = ip 3 , (19) 

c _ ■ ■ i £ 

where the scalar density is given by 

M* /- fc Fp k 2 M* /*F» k 2 

P s = ~JT L dk n.9. , ,^m/9. + ~JT L dk n.9 , »^m/s ■ ( 20 ) 



7T 2 JO (P + M* 2 )V2 ^2 y Q (£2 + M *2)l/ 

Because of the factor of 1/2 on the right-hand side of Eq. (|T§D, it follows that the relevant 
expansion parameters for the energy density are <3>/M, W/M, and 2R/M. Moreover, by 
dividing £ by M 4 and expressing the result in terms of these expansion parameters, one can 
identify the scaled couplings that should all be of roughly the same size if they are "natural" , 
namely, 

111 a ijk M i+2 i +2k -* 
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III. THE HIGH-DENSITY LIMIT 



The coupling constants c 2 , c 2 , c 2 , and in Eqs. QlSf) and ( |16D enter as unknown model 
parameters. According to the generally accepted procedure, these parameters will be chosen 
to reproduce the properties of nuclear matter near equilibrium. The basic ingredient in 
many astrophysical problems, e.g., neutron-star calculations, is the EOS 

p = p(£), (21) 

which is then extrapolated into the neutron-rich and high-density regime. Anticipating the 
results of the next section, one can expect that different parameter sets that lead to identical 
equilibrium properties produce qualitatively similar equations of state at low densities. The 
relevant question is whether this qualitatively similar behavior persists at high densities, 
particularly in the regime important for neutron stars. As a first step in this direction, 
we investigate the high-density limit of the EOS generated by the model introduced in the 
previous section. 

To make the discussion more transparent, we focus here on pure neutron matter (kp p = 
0, /cFn = kp), although nuclear matter in /3-decay equilibrium with a finite proton to neutron 
ratio is necessary to achieve accurate results for maximum neutron star masses. We will 
return to this issue in the next section. 

It is clear that a sufficiently large number of couplings introduces a high degree of flex- 
ibility. Due to the nonlinearity of the problem, not all families of parameter sets lead 
to physically acceptable results, which provides one way to restrict the parameter space. 
Classes of models can be ruled out if basic physical requirements are violated. For example, 
one certainly requires that the pressure p be a smooth function of the energy density S. 
Moreover, it is necessary that the speed of (first) sound C\ respect causality and also be real, 
to ensure stability. That is, 

In addition to these general principles, we require a positive and bounded value of the 
nucleon effective mass, i.e., 

< M* < M . (22) 

This is motivated by the expectation that physically reasonable models will demonstrate 
some degree of chiral-symmetry restoration at finite density, leading to a reduction in the 
nucleon mass, with the most extreme situation corresponding to total restoration of the 



symmetry. The consequences of Eq. (22) will become more transparent in the following. 

By specializing the formalism of the preceding section to pure neutron matter and by 
using the self-consistency equations (|l7|)-(|19l), the pressure and energy density can be ex- 
pressed as 

p = Po (k F ,M*) + Ap($,W,R) 

1 r fe F „ k 4 
ak- 



2 



3ti 2 Jo (k 2 + M* 2 ) 1 / 

+ ^ w2+ h R2 -^ 2+ ^ ~ aijk $i w2J R2k ' (23) 

v p s i,j,k 
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£ = S (k F , M*) + AS W, R) 
dkk 2 (k 2 + M* 2 ) 1 ' 2 



ii 2 Jo 

+ T2 W * + Tl R2 + $2 + E( 2 i + 2A; - !) $i ^ ^ > ( 24 ) 

^ C v ^ C p ^ C S 

where p an d £o denote the results for a relativistic, noninteracting gas of spin-1/2 baryons 
with mass M*. Note that these expressions include, as a special case, models where only 
the aojk are nonzero, so that there are no scalar-vector couplings. Moreover, the isovector 
density in Eq. (|I9"D is replaced by 

P3=~P= . (25) 

At high densities, the left-hand sides of the self-consistency equations flUf ) and (|T^) must 
grow linearly in p, and thus we start with the Ansatz 

lim W = w p a , ]imR = r ff , (26) 

p^oo p^OO 

where < a, (3 < 1, Wq > 0, and r < 0. Since we assume that the effective mass is bounded, 
we can replace the scalar field by the limit 

$00 = lim (M - M*) < M . (27) 

p— >oo 

To fulfill the resulting self-consistency equations: 

1 W + J2 2j 5«* $1 R 2k = P} (28) 

i i2 + £ 2fc a iifc ^ J**" 1 = ip 3 = -ip , (29) 

there must be integers (j m , k m ) and (j^, k' m ) with 

(2j m - l)a + 2k m (3 = 1 > (2j - 1)q + 2fc/3 for all j^j m ,k^k m , (30) 
2j> + (2*4 - l)/3 = 1 > 2ja + (2k - l)/3 for all j^j' m ,k^k' m . (31) 

(This assumes that only one term in each sum produces the leading asymptotic behavior; if 
this actually happens for more than one term in a sum, the conclusions below are unchanged.) 
Using 

Jimp = ^^p 4 / 3 + O(p 2 / 3 ) , \un£ = ^^p^ + O(p 2 ^ , (32) 
the leading contributions to the pressure and the energy density are found to be 
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lim D - ( 3?r2 ) 1/3 4/3 , J_ 2 2a , J_ 2 2/3 

p ^P- 4 P + 2c 2 oP + 2c 2 ° P 

i ,,,2j m 2fc m 1+q - , 2}' 2*4 1+/3 V" n (T)* fQQ^ 

lta £ = ^^ + ^« + ^ 

+ w?'"rj*v 1+ '»(2j„ + 2A„, -1)2 ay,„ fc „*^ 

i 

+ + 2*4 - 1) £ . (34) 

i 

To this point, the discussion is rather general. To make the conclusions more concrete, 
we discuss two distinct situations separately: 

1. No coupling between W and R. 

In this special case, the asymptotic behavior of the fields is governed by their highest 
powers in the potential @. From Eqs. (|30D and fl3jj ), we obtain 

a = — — , /3 



2 J max 1 2/c max 1 

For J max 1 , fc ma x 2^ 1 or j max > 1 , /c max = 1, the quadratic terms dominate the 
right-hand sides of Eqs. ( p3| ) and fl3~4]) , so that 

lim —W 2 oc p 2 or lim — i? 2 oc p 2 . 

This case includes the original version of the Walecka model [|J and generates the 
limiting behavior 

lim p = S . (35) 

The sums in Ap and AS contribute to the leading order only if j max = 2 , fc max > 2 
or j max > 2 , A; max = 2. In this case, the quadratic terms can be neglected, and Ap 
and AS are of the same order as the contributions from the ideal-Fermi-gas terms. 
However, the factors in Eq. (|3~4|) conspire such that 

lim Ap = - AS , (36) 

and the functional form of the limiting EOS resembles that of an ideal Fermi gas: 

lim p = -S . (37) 

p^oo 3 

In the remaining cases (j max > 2, A; max > 2), the dominant contributions arise solely 
from pq and S , which also leads to Eq. fl3"7p. Note here the importance of Eq. (P2"D, 
which implies that in the high-density limit, M* becomes negligible, at least to leading 
order. 
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2. At least one coupling between W and R. 

From Eqs fl30|) and (0), it follows directly that 

a < - and 8 < - . 

~ 3 M ~ 3 

The quadratic contributions in the fields are negligible, and Ap and AS contribute to 
the leading term only ifa = /3=l/3, where again Eq. ( jStp holds. In any event, this 
leads to the limit of Eq. 



To summarize, we conclude that the high-density limit of the EOS is strongly influenced 
by nonlinear meson-meson interactions, which agrees with the conclusion of Bodmer and 
Price PJTPfl. The limit in Eq. Q3"5| ) obtained in the original version of the Walecka model 
is a special case; in the more general situation, the nuclear matter EOS approaches that 



of an ideal Fermi gas, given by Eq. (|37f) . We will show in the next section that these two 
limits can be achieved using different models with parameter sets that reproduce the same 
equilibrium properties of nuclear matter. 



IV. CONSEQUENCES FOR NEUTRON STARS 

The high-density limit of the EOS and the way in which the matter approaches the 
asymptotic regime have important consequences in neutron star calculations. The masses 
and radii of stars are sensitive to the stiffness of the EOS, thus providing a quantitative 
measure for studying the impact of the nonlinear interaction terms in Eq. 

To be specific, it is necessary to choose an explicit potential, and we begin with the form 

E CW ( V^) fc = - + ^giW) 3 + Vb") a , (38) 

i,j,k ■>. . . 

which includes a subset of the meson self-interactions up to fourth order in the fields. As 
discussed in the Introduction, setting some of the allowed cubic and quartic couplings to 
zero is "unnatural" , but as we will discover, the model defined by Eq. (|38|) is already general 
enough to produce significant differences in predicted neutron star masses, and restoring the 
omitted couplings will lead to even more variation in the results. Moreover, the present model 
can be related to the most common models discussed in the literature, and it generalizes 
them to include a nonlinear isovector interaction. The motivation for adding the quartic 
rho-meson term is that one expects this coupling to be essentially unconstrained by normal 
nuclear observables, where the neutron-proton asymmetry is low, but it may have significant 
impact on the neutron-rich matter in neutron stars. As noted earlier, since the meson fields 
are effective (Kohn-Sham) potentials, we are not concerned here with their elementary 
excitations, and considerations of causality [0,0 are unimportant. 



In nuclear matter calculations, this model contains seven free parameters. The poly- 
nomial in Eq. fl38|) contains four couplings that we may write as R = R/gl, A = X/g^, C> 
and £; in addition, values for the three ratios c 2 = gf/mf (i = s, v, p) are needed. Five of 
the seven parameters can be chosen to reproduce the equilibrium properties of symmetric 
nuclear matter, which we take as the equilibrium density and binding energy (p , — eo), the 
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nucleoli effective (or Dirac) mass at equilibrium (Mg), the compression modulus (Kq), and 
the bulk symmetry energy (a^). The first three of these are tightly constrained ||11|| , whereas 
the latter two are not. In principle, the sensitivity of the high-density EOS to reasonable 
variations in K Q and 04 could be examined, but for simplicity, we keep their values fixed in 
most of our calculations. The "standard" set of equilibrium properties used here are listed 
in Table [IJ these are motivated by successful descriptions of bulk and single-particle nuclear 
properties fTTJil,2~I|. The nucleon mass is fixed at its empirical value (M = 939 MeV). 



Our primary goal is to study the influence of the nonlinear vector-meson interactions 
on neutron star masses. Since there are more free couplings than normalization conditions, 
we proceed as follows: We choose values for the couplings ( and £ and determine the 
remaining couplings by requiring that they reproduce the desired equilibrium properties. 
This is achieved by solving a set of transcendental equations that relate the parameters 



directly to the nuclear matter properties |19| , 21| . Although we have no specific guidance 



on the allowed values of ( and £, we rely on the assumption of naturalness, and based on 
the discussion at the end of Sec. II, we observe that the following parameter combinations 
should all be of roughly equal size: 

1 1 1 R \ ( f 

and 



2c 2 M 2 ' 2c 2 w M 2 ' 8c 2 p M 2 ' 6M ' 24 ' 24 ' 384 ' 

Typical values for the first three parameters are between 0.001 and 0.002, so that the natural 
values of ( and £ are roughly limited to < ( < 0.06 and < £ < 1.0. (To avoid abnormal 
solutions of the vector field equations, i.e., those with finite mean fields at zero density, £ 
and £ must be positive. 4 ) We will include results for vanishing ( and £, which are in a strict 
sense unnatural, in order to make contact with earlier calculations. For a more thorough 
discussion of naturalness and its implications, see Ref. |^T] . 



Using the notation of Sec. II, the self-consistency equations (|T^)-(|i~9"|) can be written as 

1 $ + l& + V = Ps , (39) 



c 2 2 6 




(40) 



2 P3 = ~2 P ■ (41) 



The expressions for the pressure and the energy density follow as 

1 /•** , r k 4 
P = 7; — 9 / dk- 



3n 2 Jo (k 2 + M* 2 y/ 2 

+ -L W i + 1 W * + * R 2 + 1 B 4 - -At $ 2 - -$ 3 - — $ 4 , (42) 
2c 2 24 2c 2 p 24 2c 2 6 24 ' v ; 



4 Note that positive £ and £ imply that the resulting nonlinear interactions are attractive. This 
constraint on the highest-order vector interactions appears to be general and explains why our 
earlier analysis finds that nonlinear interactions soften the equation of state. 
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£ 1 [ kF dkk 2 (k 2 + M* 2 )^ 2 



+ ^H/2 + C^4 + 1 i?2 + i jR 4 + 1 $2 + ^ $3 + A$4 . (43) 

2 C 2 8 2c 2 p 8 2c s 2 6 24 v 7 

We begin our discussion with the model introduced by Bodmer and Price 0, which 
corresponds to £ = 0. According to the discussion in the preceding section, this model has the 
interesting feature that the high-density EOS of pure neutron matter approaches p = £, while 
in symmetric matter, where the mean-field of the p meson vanishes, the EOS approaches 
the massless Fermi gas limit, given by Eq. (|37|) . In Fig. [I], we show the binding energy 
curves for symmetric and pure neutron matter for different values of the nonlinear coupling 
£. We emphasize that all parametrizations reproduce the same equilibrium properties listed 
in Table |. 

At low densities, all the curves approach a common limit, because the nonlinear terms 
do not contribute at leading order in a low-density expansion. At higher densities, the 
softening of the EOS as ( increases is clearly visible, at least for symmetric matter. The 
softening in neutron matter is more apparent in Fig. |^. In the regime of intermediate density, 
200 < £ < 1000 MeV/fm 3 , the EOS becomes softer with increasing values of ( |19|]. To study 



the approach to the asymptotic limit, we examine the nonleading terms in the high-density 
expansion: 

2/3 

lim p = £ - ( - ' - ) £ 2 ^ + 0(£ 1 / 3 ) for C = 0, (44) 



£— >oo 



lim p = £ 




1 



£ 2/3 + 0(£ 1/3 ) for Ct^O, (45) 



£-mx> r Qtt 2 

which reveals two important features. First, the C ^ results are nonanalytic in (. Thus 



one cannot reproduce Eq. (fH) by taking the £ — > limit of Eq. fl4"5j ); at least as far as £ 
is concerned, the high-density expansion is essentially a strong coupling expansion. Second, 
we observe that the coefficient of the nonleading term is smaller for ( = for two reasons: 
the appearance of the isoscalar coupling c 2 in the denominator and the absence of the 
multiplicative factor containing (. (Note that c 2 is independent of (.) These two features 
produce a coefficient that is roughly an order of magnitude smaller for ( = than for ( ^ 0, 
which explains the relatively slow approach to the asymptotic limit in the latter case, as is 
evident from Fig. 0. 

The consequences for neutron stars can be studied in Fig. |||, where the star masses are 
shown as a function of the central mass density p c . As expected from Fig. |2|, the maximum 
mass decreases with increasing £. This decrease is substantial: from M max = 2.9 M & for 
( = to M max = 2.1 Mq for ( = 0.06, which is roughly 30%. The shifts in the maximum 
mass are most dramatic for small couplings; for larger couplings, the softening effects begin 
to saturate. 

To understand this result, it is useful to identify the regime of energy density that is most 
important in determining the mass of the star. This regime can be deduced from Fig. |], 
where we show the radial mass density distributions for several neutron stars, as well as the 
corresponding energy densities. Observe that most of the mass is generated at radii between 
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6 and 12 km, which corresponds to energy densities of several hundred MeV/fm 3 . As can 
be seen from Fig. ^], this includes the regime where the EOS is sensitive to (, because this 
is where the W A contribution to the EOS begins to become important. Note, however, that 
the contribution of the W A term in this regime is still smaller than that of the W 2 term; the 
quartic term does not begin to dominate until the energy density reaches several thousand 
MeV/fm 3 , as indicated by the coalescence of the dashed and dotted curves in Fig. ||. 

We now return to the general form in Eq. ( |3"8| ) and allow nonlinear self-interactions of 
the p mesons. In contrast to the isoscalar coupling £, the new quantity £ does not enter 
in the calculation of our five "standard" equilibrium properties of nuclear matter, and the 
other parameters are determined independently. This is true even for the symmetry energy, 
because the new coupling £ first appears at order (N— Z) A in an expansion around symmetric 
nuclear matter. In principle, £ could be constrained by fits to liquid-drop expansions of the 
energy, but in most such fits this parameter is set to zero. (We found only one nonzero 
value in the literature pSfl .) Thus, at present, contributions to the symmetry energy beyond 
terms of order (N — Z) 2 are practically unconstrained. 

Fig. |5| shows the binding energy of nuclear matter as a function of the proton fraction y 
for two different densities. The curves are calculated for various values of £ with ( held fixed. 
For clarity, we plot the fractional shift in S relative to its value with £ = 0. At normal nuclear 
density, results for different £ are virtually indistinguishable, but at high density, the curves 
differ by a few percent when the proton fraction becomes very small. This demonstrates that 
it is possible to generate families of models that reproduce identical properties of nuclear 
matter at low and normal densities, but which generate different predictions at high densities. 

According to the analysis in the preceding section, the EOS asymptotically approaches 
the massless Fermi gas limit for both symmetric and neutron matter. The high-density 
expansion for neutron matter corresponding to Eq. is 



1 n 2vr 
hm p = —c H 

S^co 3 3 



W_2_\ 2 / 3 , J_/J_N 2 / 3 



i + m 1/3 + (^ 1/3nV2 



£i/2 + o(£°) . (46) 



In contrast to the previous case, the asymptotic limit is now approached from above. More- 
over, the analytical form has changed, since Eq. (EH) indicates a series in powers of E 1 ! 2 . 

The dependence of the neutron star masses on £ follows the same trend as found for £. 
This can be gleaned from Fig. [(], which shows results for pure neutron matter. For fixed (, 
the maximum mass decreases with increasing £, although the effect of this new parameter is 
smaller (less than a 10% change in the maximum mass), given the limitation on parameter 
values imposed by naturalness. 

A more complete picture of M max is given in Fig. ff], where the variations with both ( 
and £ are shown. It is apparent that parameter sets that yield identical properties near 
nuclear equilibrium can still generate values of the maximum mass that differ by as much 
as one solar mass. We find variations between M max = 2.9 M and M max = 1.9 M Q for stars 
composed of pure neutron matter. The masses are most sensitive to changes at small values 
of the couplings, particularly for £ m 0, which can be related to the nonanalytic form of the 
EOS in terms of ( and £. [See Eqs. (fH])-(|46|).] As seen earlier, the effect of £ is smaller than 
that of (. 
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So far, we have considered pure neutron matter, which gives only a qualitative picture 
of neutron star properties. For a more realistic description, it is necessary to consider beta- 
stable matter, as this includes the protons, which softens the EOS. Maximum neutron star 
masses for beta-stable matter are also shown in Fig. |7]. The dependence on the isoscalar 
coupling ( is similar to that obtained earlier (the maximum mass varies between M max = 
2.8 M Q and M max = 1.8 M ), but the influence of the isovector coupling £ is less drastic, 
since the matter becomes significantly more symmetric in the region that gives the largest 
contribution to the mass. 

In Fig. |^, we examine the dependence of the maximum mass on the compression modulus 
at equilibrium, which is not particularly well known. (Relativistic mean-field models with 
200 ^$ Kq < 350 MeV can produce accurate nuclear binding-energy systematics and surface 
energetics pT| .) The shaded band shows the total predicted variation in maximum mass 
when both ( and £ are varied within the bounds imposed by naturalness. The dashed curve 
shows the predicted variation when K Q is varied at fixed ( and £. Evidently, the variations 
in the maximum mass arising from the possibility of nonlinear vector meson interactions is 
much greater than that arising from the uncertainty in the compressibility. 

These results raise the interesting question of whether the only significant nonlinearity 
is the quartic, isoscalar vector interaction. In other words, once one takes £ ^ to soften 
the EOS at high densities, does the addition of further nonlinearities produce only small 
effects? To examine this question, we extend the model of Eq. (|38|) to include sixth-order 
and eighth-order terms involving the vector-isoscalar meson: 

AV = - C -/AV^f - C ^g 8 v (V,vn 4 ■ (47) 

Thus C'/6! and ("/8\ are the relevant ratios to be included in the parameter list given earlier, 
and we will initially set £ to zero and examine the consequences of varying (, and (" 
within the bounds imposed by naturalness. 

Figure ^ shows neutron star masses as higher-order vector nonlinearities are included 
sequentially. Evidently, the quartic interactions are the most important, producing a roughly 
30% variation in the maximum. The effects of the sixth-order term are quite modest (roughly 
10%), while the eighth-order contributions are essentially negligible (roughly 2%). Here the 
parameters are varied within the natural ranges < £ < 0.06, < £' < 1.2, and < Q" < 60. 
Thus we have the encouraging result that once the W A interaction has been accurately 
calibrated, contributions from higher-order nonlinearities are relatively unimportant. To 
indicate the most extreme reduction in maximum mass possible in the present model, the 
cross in Fig. |9| shows M max = 1.58 M , which is obtained for beta-stable matter when the 
couplings ( = 0.06, (' = 1.2, (" = 60, and £ = 1.5 are included. Note that this value of 
M max is only slightly (~ 10%) larger than that of the most massive observed neutron stars. 

Our analysis to this point has revealed significant model and parameter dependence in the 
high-density EOS. It is therefore of interest to see if these variations influence predictions 
arising from other relevant dynamics in systems with high densities. As an example, we 
study the effect of the high-density hadronic EOS on the existence of quark-matter cores 
in neutron stars. We adopt a simple two-phase model |37],[]J based on a first-order (van 
der Waals) phase transition between the hadronic and quark phase. Although there are 
indications from QCD lattice calculations that the hadron/quark phase transition is second- 
order at vanishing chemical potential [0,0] (for massless quarks), the true behavior of the 
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transition at finite density (and indeed, whether one actually exists) is unknown at present. 
Until more reliable information is available, one must resort to a separate description of the 
quark phase and the hadronic phase. Moreover, whereas a detailed description involving 
beta-stable matter requires a careful treatment of the phase transition in systems with two 
conserved charges (baryon number and isospin) |27|,[E|, here we will be satisfied with a 
qualitative discussion based on pure neutron matter. This is certainly reasonable, given the 
large uncertainties we have already found in the hadronic EOS at high density. 
We adopt the simple EOS involving massless u and d quarks given by 

V = \s - \b , (48) 

where the confinement property of QCD (or alternatively, the anomaly in the trace of the 
energy-momentum tensor) is modeled by a positive constant b, which represents the energy 
per unit volume in the vacuum. 

We return to the hadronic model of Eq. (|38|), which leads to the high-density expansion 
in Eq. (f46|). Our discussion is based on the simple observation that independent of the actual 
nature of the transition, it is driven purely by the energetics in the two phases. 

By comparing Eq. fl4*8|) with the expansion in Eq. (f46|) , one observes that in the quark 
phase, the limit p = S/3 is approached from below, whereas in neutron matter, for the 
general case ( > and £ > 0, the limit is approached from above. To decide whether a 
transition takes place, it is necessary to compare the energy/baryon in both phases. For the 



quark phase it is expressed as [37 



£/p=\^f{a s )p^ + b/p, (49) 

with 

/(a s ) = (l + 2 4 / 3 )(l + ^), 

which includes the lowest-order contribution in a s (the exchange energy). As discussed in 
the preceding section, in hadronic models that are characterized by Eqs. ([44]) and (0), the 
quadratic terms dominate at high densities, so that 

lim £ I p (x p , (50) 

for C > and £ = 0. This is also true if ( = and £ > 0. Thus, at sufficiently high 
densities, neutron matter always has higher energy compared to the quark phase P? |, and 



the two phases can be connected by a Maxwell construction, which signals the transition 
from hadron to quark matter. 

The situation is different in the general case ( > 0, £ > 0. For the asymptotic form of 
the energy corresponding to Eq. (flop, one obtains 



lim£/p=|(37rY /3 



/ 2 n 1/3 1 / 1 \ 1/3 

1 + fec) 



?l* + 0{p-V>), (51) 



which is, up to the prefactor, the same leading behavior as in Eq. (|49|) for the quark phase. 
Therefore a phase transition is possible, i.e., neutron matter has higher energy, only if 
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3 l/3 



1 + 



2x1/3 l/lxl/3 



- 



> f(a B ) ■ (52) 



This remarkable observation implies that for sufficiently large values of the nonlinear cou- 
plings, the matter remains in the hadron (neutron) phase, at least in the simple model 
discussed here. More generally, one observes that independent of the asymptotic form, in- 
creasing the couplings £ and £ increases the density at which the phase transition occurs (if 
it does), since the hadronic EOS becomes softer as the nonlinear couplings increase. 



The different possibilities are illustrated in Fig. [10]. The curve labeled a corresponds 
to Eq. (|50D ; the asymptotic behavior is clearly different from the quark EOS. The phase 
transition occurs at roughly 3po, and the Maxwell construction is indicated by the dotted 
line. In the situation described by curve b, neutron matter and quark matter have a similar 
asymptotic behavior, but the condition Qo*^ ) remains true, and the two curves cross, leading 
to a phase transition at roughly 6p . Finally, curve c lies completely below the quark EOS 
and a transition is not possible; neutron matter is stable at all densities. 

If one introduces interaction terms of higher than fourth order in the fields, for example, 
{V^V^) 3 , the hadronic energy is dominated by the Fermi gas contribution, and Eq. (|5lD must 
be replaced by 

lim£/p= ^STr 2 ) 1 ^ + (p-W). (53) 

Since 3 1 / 3 < f(a s ), no transition is possible in the asymptotic regime for any choice of 
hadronic parameters in this case. On the other hand, the transition regime also depends 
on the model and parameters used for the description of the quark phase. In our model, 
the vacuum constant b and strong coupling a s can be used to shift the transition point 



substantially ||52| , p3| , so that a transition may occur outside the asymptotic regime. It is 
clear, however, that the strong model dependence in the hadronic EOS introduces large 
uncertainties in any attempted prediction of these values. 

The consequences for neutron stars in this model follow straightforwardly. Nonlinear 
vector meson interactions soften the hadronic EOS, which lowers the maximum neutron star 
mass and increases the density of the transition to quark matter. In contrast, a stiff hadronic 
EOS lowers the density of the phase transition, and since the quark matter EOS is soft, also 
tends to decrease the maximum star mass. Thus it may be impossible to decide, from 
neutron star masses alone, whether quark matter cores exist in neutron stars, and similar 
conclusions may be drawn about other exotic phenomena that soften the EOS. (The situation 
is complicated further by the continuous nature of the transition when two conserved charges 
are involved, which is the more physical case [0.) Whereas it might be possible, using the 
results of more advanced calculations of the finite-density hadron/quark phase transition, 
to rule out certain parametrizations of the hadronic EOS, existing uncertainties in both the 
nature of the phase transition and in the high-density hadronic EOS preclude any definite 
conclusions at this time. 



V. SUMMARY 

In this paper we study the equation of state of nuclear and neutron-star matter based on 
relativistic mean-field theory. Our starting point is an effective action (or energy functional) 
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containing Dirac nucleons and local scalar and vector fields. These fields are interpreted 
as relativistic Kohn-Sham potentials, and nonlinear interactions between the fields are in- 
troduced to parametrize the density dependence of the energy functional. We calibrate 
the energy functional by observing that at normal nuclear densities, the ratios of the mean 
fields to the nucleon mass are small, and thus the nonlinear interactions can be truncated 
at some low order in the fields. The unknown parameters can then be fit to properties of 
nuclear matter near equilibrium that are known to be characteristic of the observed bulk 
and single-particle properties of nuclei. 

We then extrapolate the resulting equation of state into the neutron-rich, high-density 
regime to calculate the properties of neutron stars. Two problems arise in the extrapolation: 
First, even with a significant truncation of the energy functional, the unknown parameters 
are underdetermined. Thus there exist families of parameters that reproduce exactly the 
same nuclear matter properties near equilibrium, but which produce potentially different 
high-density equations of state. Second, terms omitted from the functional because they are 
negligible at normal density may become important at densities relevant for neutron stars. 
This is true even if we assume that the coupling parameters are "natural" , which means that 
they are all of roughly the same size when expressed in appropriate dimensionless ratios. 

Our basic goal is to determine, in light of these two problems, whether the calibration 
at equilibrium nuclear matter density is sufficient to predict a maximum neutron star mass 
within a reasonably small range. This is relevant in view of recent calculations that hope 
to see evidence for "new" physics in neutron stars (such as quark cores, strange matter, or 
kaon condensates) based on the need for a softer high-density equation of state than that 
provided by neutrons, protons, and electrons alone. These calculations assume that the high- 
density behavior of these more mundane components is well known, and in particular, that 
the contributions of many-nucleon forces are negligible [27-^. These many-body, density- 
dependent forces are precisely the ones that are difficult to calibrate using observed nuclear 
properties; the question is whether one can build nuclear equations of state with different 
types of many-body forces that all reproduce the observed properties near equilibrium, but 
which yield significantly different results at high density. 

By beginning with a meson self-interaction potential containing arbitrary powers of 
scalar-isoscalar (a), vector-isoscalar (uj), and vector-isovector (p) fields, we show that the 
meson nonlinearities can have a profound effect on the high-density equation of state. In 
models where the vector mesons enter the potential at most quadratically, the equation of 
state is stiff and asymptotically approaches p = £. (The Walecka model is a special case.) 
In models where the vector fields enter with high powers, these fields become negligible at 
high density, and the asymptotic equation of state resembles that of a free, relativistic gas: 
p = S/3. The intermediate case occurs when the vector fields enter quartically; the asymp- 
totic equation of state is still soft (p = S/3), but the approach to the asymptotic limit is 
determined by the coupling parameters. 

We illustrate these results using specific models containing quartic uj and p meson cou- 
plings and also sixth- and eighth-order uj couplings. All models are calibrated to exactly the 
same nuclear properties at equilibrium, for all choices of parameters. We find that by far 
the most important coupling constant is that of the quartic uj term; even when restricted 
by the requirements of naturalness, variations in this parameter can produce variations of 
nearly one solar mass in the predicted maximum neutron star mass. (This is true for both 



19 



pure neutron matter and for beta-stable matter.) This uncertainty is clearly relevant on the 
scale on which one hopes to identify new, exotic effects. Moreover, increasing the strength of 
the quartic coupling softens the equation of state, which is precisely the effect sought from 
the exotica. The impact of a quartic p meson coupling is smaller, and its effects are only 
appreciable in stars made of pure neutron matter; the maximum masses of stars computed 
with beta-stable matter show little change when this parameter is varied within the bounds 
imposed by naturalness. Similarly, sixth-order uj interactions have only a modest effect on 
the predicted maximum mass, and by the time the eighth-order terms are included, the 
high-density equation of state is already so soft that these terms are negligible. 

We emphasize that the importance of these many-body effects is not limited to the 
domain of relativistic mean-field theories; equations of state based on nonrelativistic poten- 
tials P6| use interactions that are calibrated to few-body systems and that are insensitive 
to possible six- or eight-body forces that may be relevant at high density. This is especially 
important because a mean-field calculation with four-component spinors and just two-body 
Lorentz scalar forces already implicitly contains an infinite string of many-body forces if it 
is recast in terms of two-component spinors. 

To illustrate these difficulties more concretely, we also study the role of the hadron/quark 
phase transition in a simple two-phase model with a first-order transition. We find that with 
couplings well within the bounds of naturalness, it is possible to push the phase transition 
to arbitrarily high density, and even to make it disappear altogether. Although the absence 
of a phase transition is probably unrealistic and could serve to exclude some values of the 
hadronic parameters, there are still too many uncertainties on both sides of the transition 
(as well as in the nature of the transition itself!) to make any definitive statements. 

We therefore conclude that existing methods for calibrating the nuclear equation of state 
for extrapolation into the neutron-rich, high-density regime appropriate for neutron stars 
cannot constrain the predicted maximum star mass well enough to make reliable statements 
about the existence of "new" physics beyond the dynamics of neutrons, protons, and elec- 
trons. We show that the uncertainties arising from an incomplete knowledge of the hadronic 
many-body forces are much larger than those arising from an imperfect knowledge of the 
properties of nuclear matter near equilibrium, such as the nuclear matter compression mod- 
ulus. Even the old question of the role of the hadron/quark phase transition is problematic, 
since the high-density hadronic equation of state can be made essentially as soft as desired 
by the addition of nonlinear interactions that are still consistent with equilibrium nuclear 
matter properties. 

One positive conclusion is that the most important nonlinear parameter is that of the 
quartic u interaction. If this term could be accurately calibrated, the uncertainties intro- 
duced by other nonlinear interactions are likely to be tolerable. (One caveat: we do not 
study carefully the sensitivity to variations in mixed scalar-vector interactions because of 
the overwhelming sensitivity to the quartic uj term; if the latter were well constrained, the 
role of scalar-vector interactions should be examined in more detail.) Although this coupling 
has not been extensively studied in mean-field calculations, there are several possibilities for 
determining it reasonably well. First, since a quartic uj interaction leads to a nonlinear 
density dependence in the vector part of the baryon self-energy, one could calibrate this 
interaction by fitting to the self-energy obtained in a Dirac-Brueckner-Hartree-Fock calcu- 
lation, for example. Although some initial attempts at this procedure have been made PO 
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the resulting parameters are not always natural; it is probably necessary to repeat the pro- 
cedure using all possible scalar and vector self-interactions through fourth order and to fit 
both the scalar and vector part of the self-energy simultaneously. Second, the nonlinearities 
in the energy functional can be interpreted in terms of effective masses for the vector and 
scalar mesons (defined by diagonalizing the matrix of appropriate second derivatives of the 
energy functional). This may provide useful constraints in the future, if concrete empirical 
information on these effective masses becomes available. Third, it is possible that additional 
observables in finite nuclei could constrain the nonlinear interactions. For example, some 
recent work suggests that the ratio of the nuclear matter "skewness" (which is related to 
the third derivative of the energy with respect to density at equilibrium) to the compres- 
sion modulus K is constrained by monopole vibrations. Although a recent calculation of 
nuclear ground-state properties shows little correlation with this ratio [^TJ, a more detailed 
examination of dynamical effects could provide meaningful constraints. 

To summarize, precise predictions of the properties of neutron stars apparently require 
more accurate calibrations of the nuclear equation of state than are currently available. It 
is especially important to have the high-density behavior of the "standard" components 
(neutrons, protons, and electrons) under control before one can make reliable statements 
about the existence of "new" physics. Since the window on experimentally observable nuclear 
properties is a narrow one, producing an equation of state that can be extrapolated with 
confidence remains a major challenge. 
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FIGURE CAPTIONS 



FIG. 1. Binding energy of symmetric and neutron matter. 

FIG. 2. Equation of state of neutron matter. 

FIG. 3. Neutron star masses for different nonlinear couplings (. Pure neutron matter is 
assumed. 

FIG. 4. Density distributions of neutron stars as a function of the radius, calculated with 
different parameters £. The central mass density p c = 1.5 x 10 15 g/cm 3 is the same for all curves. 
Part (a) shows the radial mass densities and part (b) the corresponding energy densities. The stars 
acquire most of their mass from the region 100 < £ < 500 MeV/fm 3 . 

FIG. 5. Energy difference A£ = [£(£ = 0) — £(£)]/£(£ = 0) of nuclear matter as a function of 
the proton fraction y calculated at constant baryon density. In part (a) the density is fixed at its 
equilibrium value po and part (b) shows the results at 5po- Note the different vertical scales in (a) 
and (b). 

FIG. 6. Neutron star mass for different values of £ at fixed (. Pure neutron matter is assumed. 

FIG. 7. Maximum neutron star mass as function of £ and £. Results for pure neutron matter 
and for matter in /3-equilibrium are displayed. The shaded areas show the mass range obtained 
when £ is varied; the upper boundaries correspond to £ = and the lower boundaries to £ = 1.5. 

FIG. 8. Maximum neutron star mass as a function of the compression modulus Kq. (All 
other nuclear matter inputs are held fixed.) The shaded area marks the covered range of masses. 
The upper boundary corresponds to £ = 0, £ = and the lower boundary to £ = 0.06, £ = 1.5. 
For fixed values of the nonlinear couplings, the mass changes marginally with the compression 
modulus. This can be seen from the dashed line inside the shaded area, which corresponds to 
C = 0.02, £ = 0.5 

FIG. 9. Neutron star masses for models with different nonlinear couplings for the neutral 
vector meson. The uppermost curve corresponds to the Walecka model, including cubic and quartic 
couplings for the scalar meson. Region A shows the range of masses obtained when the quartic 
vector coupling is turned on. Regions A and B correspond to models with up to sixth-order terms, 
and regions A, B, and C include an eighth-order term. The values of the nonlinear couplings 
are chosen within a natural range, as described in the text. The cross indicates the maximum 
mass (1.58 M ) obtained for beta-stable matter in a calculation that includes the £, £", and £ 
couplings (see the text). 
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Energy per baryon for neutron matter (solid) and quark matter (dashed) . The quark 

120 MeV/fm 3 and a s = 0.4 



FIG. 10 

matter results are calculated using b 
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TABLES 

TABLE I. Equilibrium Properties of Nuclear Matter 



(k F )o 


Po 


M */M 


e 


K 


04 


1.30 fm" 1 


0.1484 fm- 3 


0.60 


-15.75 MeV 


250 MeV 


35 MeV 
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